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We propose a new formalism to handle asymmetric beams in the data analysis 
of cosmic microwave background anisotropy experiments. For any beam shape, 
the formalism finds the optimal circularly symmetric equivalent and is thus eas- 
ily adaptable to existing data analysis methods. We demonstrate certain key 
CN ! points by using a simulated highly elliptic beam, and the beams and data of the 

P^ ! MAXIMA- 1 experiment, where the asymmetry is mild. We show that in both 

^ ! cases the formalism does not bias the angular power spectrum estimates. We 

O ' analyze the limitations of the formalism and find that it is well suited for most 

^—^ ' practical situations. 

Q . Subject headings: cosmic microwave background — cosmology:theory — large-scale 

^ , structure of the universe — methods:numerical 
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INTRODUCTION 



A new generation of Cosmic Microwave Background (CMB) mapping experiments is 
beginning to produce data of unprecedented quality (see e.g., Torbet et al. 1999; Miller 
et al. 1999; DeBernardis et al. 2000; Hanany et al. 2000, hereafter HOO). Much of the 
experimental effort is concentrated on probing angular scales of about 10 arcminutes. To 
fully benefit from the scientific potential of these high resolution data sets, a new level of 
sophistication is required in quantifying all possible sources of error in the experimental 
procedure and data analysis pipeline (e.g., Ferreira & Jaffe 2000). Particular care must 
be used to accurately quantify the instrument response to the signal and to include such 
response in the data analysis. 

In all analyses of CMB data so far the experimental beam has been assumed to have a 
radial symmetry. This assumption has been incorporated in most map-making and angular 
power spectrum {Cf) estimation algorithms (e.g.. Bond, Jaffe, & Knox 1998) and is necessary 
because of limitations in computing capability. A crude symmetric-beam approximation was 
adequate in the past since most of the error budget was dominated by statistical and other 
systematic uncertainties. However, with the precision of current and future analyses, it 
becomes essential to establish a methodology for accurately quantifying the degree of beam- 
asymmetry and properly incorporating it into the data analysis pipeline. If the beam is 
incorrectly incorporated in the data analysis pipeline, one may not only artificially distort 
the underlying structure of the measured CMB signal but also bias the estimate of the CMB 
angular power spectrum. In this paper we present a new formalism for estimating the power 
spectrum that can handle any beam shape. We show that the formalism can be applied to 
a broad variety of cases which encompass most practical applications. As a consequence, 
the detailed shape of the antenna beam should no longer pose a limitation in measuring the 
angular power spectrum of CMB experiments. 

The asymmetry of beams may arise from a variety of sources. For example, it may 
be due to the optics, or due to the finite response time of a detector which leaves imprints 
in the direction of the scan (e.g., Hanany, Jaffe, & Scannapieco 1998). Regardless of the 
origins of the asymmetry, the framework we shall present is general, and consists of finding 
an equivalent symmetric beam that replaces the asymmetric beam in the analysis of the 
data. 

Using the formalism one can assess the degree of asymmetry of a beam (see eq. [3-4]), 
how the asymmetry propagates through the analysis pipeline, and how to find an azimuthally 
symmetrized beam that best approximates the asymmetric beam (see e.g., eq. [5-9]). The 
symmetrized beam is then used in the symmetric-beam approximation of the C^ estimation 
(see eq. [5-11]). The formalism quantifies the errors introduced in the C^ estimates because 
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of the use of the symmetrized-beam approximation, the uncertainty in the final Ci estimates 
resulting from the uncertainty in the beam measurement (see eq. [7-5]), and the smoothing 
effects due to the pixelization of the map (see eqs. [8-8] and [8-9]). It also shows how to 
combine beams from independent experimental photometers (see eqs. [3-5], [4-10], [4-11], 
and [4-12]). Some useful conditions under which this new formalism will be needed are also 
provided (see eqs. [6-5] and [6-7]). 

The structure of this paper is as follows. In section 2, we describe the framework of 
CMB data analysis for the estimation of the power spectrum, so as to illustrate the problems 
related to asymmetric beams. In section 3, we define the 'index of asymmetry' (lOA) wi, 
a useful parameter in quantifying the level of asymmetry of a beam. Similarly, we define 
the 'index of combined asymmetry' (lOCA) W^, which is useful when combining data from 
photometers of different beam shapes. In sections 4 and 5, we investigate the problems 
associated with asymmetric beams. We introduce the 'average pixel-beam expansion', Bp£„i, 
and the 'pixel-pixel beam expansion', Bl^g., to provide an approximation scheme where the 
convolution effect of asymmetric beams is treated as circularly symmetric. The biasing effects 
of this approximation in the resulting estimated power spectrum Ce are also considered. In 
section 6, we derive the conditions under which one needs to employ the new formalism for 
treating asymmetric beams. In section 7, we investigate the uncertainties in the Ci estimates 
resulting from the uncertainties in the measurement of beam shape. In section 8, we discuss 
another convolution effect due to the pixelization of the CMB map. Although this is not 
a beam-related issue, we demonstrate a simple way to incorporate its treatment into our 
framework. In section 9, we numerically verify certain key points developed in sections 3 to 
8, as well as the accuracy of the proposed approximation in treating asymmetric beams. In 
particular, we use the data from the MAXIMA-1 experiment as an example to demonstrate 
the generic treatment of asymmetric beams in CMB experiments. It is shown that our 
formalism has no biasing effects in the resulting Ce estimates. Finally in section 10, we 
summarize the procedure in applying our formalism to experiments, discuss its availability, 
and draw a conclusion. 



2. THE CONVENTION AND PROBLEMS 

We first consider the standard procedure for the power spectrum estimation. This 
consists of two main steps. First, one estimates the pixelized map rUp from a given time- 
stream dt, i.e., to translate the observation from the temporal (t) to the spatial {p) domain. 
Second, one estimates the power spectrum Ce from the map nip. 
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In the temporal domain, what we observe is 

dt = -ft + nt, (2-1) 

where 74 is the CMB signal and rit is the instrumental noise. Traditionally we model the 
CMB observation as 

7t = AtpSp, (2-2) 

where we use the Einstein summation convention here and below when appropriate (usually 
over pixels and time samples, but not over spherical harmonic indices). Here Atp is the 
pointing matrix giving the weight of pixel p in observation t, and Sp is the CMB signal on 
the pixel convolved by a pixel beam -Bp(x): 

00 e 

l=Q m=-e 

where Y^^ are the spherical harmonics, and Bp£m and a^^ are the multipole expansions of 
Bp(x.) and the CMB signal respectively. Note that we use a two-dimensional vector x to 
denote locations on the surface of the sphere, which we shall often consider in the small-field 
limit (see later). 

We usually take the pointing operator Afp to be one when observing pixel p at time 
t and zero otherwise. That is, we model the signal 7^ to be the same for any observation 
within pixel p. In effect, we take the sky to be smoothed with a top-hat of shape given by 
the pixel boundary. We shall see in section 8 that, as expected, this is equivalent to an extra 
convolution included in Bp. 

With this modeling, one can thus estimate the pixelized map from the temporal data. 
This involves maximizing the likelihood of the signal given the data: 

£(s) oc Prob[rf|s] = {2ti)'^'/^ x 
expl--{n^N'^n + Tr[\nN])\ , (2-4) 

where d = dt, s = Sp, and n = rit, all as defined in equations (2-1) and (2-2), A/^ is the size of 
the time-stream, and A^ = Nu' = (ntuj,) is the time-time noise correlation matrix. Here we 
have assumed that the noise is Gaussian and that all CMB maps are a priori equally likely. 
Maximizing over s gives 

mp = m = {A'^N-^A)-\A^N-^d) 

= Sp + Up, (2-5) 
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where A = Atp as defined in equation (2-2) and np is the noise in the pixel domain. 

One then moves on to estimate the power spectrum of the map, C^ = (|a^mP)- This 
requires the maximization of the hkehhood function 

C = PToh[m\Ce\ = (27r)--^'=^/2 ^ 

exp i — [m^M'^rn + Tr [In M]) I , (2-6) 

where Mc^ is the dimension of the parameter space of C^, and 

M = Mppi = Cspp' + Cnpp', (2-7) 

with 



Cspp' — 2_^ Z-^ CeBpimB*,(. 



oo 

e=o m=-e 

F,„(xp)F;„(xpO, (2-8) 

Cnpp' = {A^N''A)-\ (2-9) 

Here Cspp' = {spSp,) is the pixel-pixel CMB signal correlation matrix, and Cnpp' = {npUp,) 
is the pixel-pixel noise correlation matrix. 

We note first that in the estimation of C^, although there exists methods like the 
quadratic estimator (Bond et al. 1998) which avoid a direct evaluation of equation (2-6), 
the relationship between the beam expansion Bpim and the power spectrum Ce remains the 
same and is illustrated in equation (2-8). Second, if the beam is identical for all pixels and 
circularly symmetric, i.e., Bpim = Bpr^j^ = Bi, then equation (2-8) can be greatly simplified 

as 

v^ 2i + 1 
Cspp' = J2 -^CiB^iPiicos Opp,), (2-10) 

£=0 ^ 

where Pi is the Legendre function and 6ppi = |xp — Xp/| is the angular distance between the 
pixels. 

Generally it is impractical to estimate C^ for all i due to the constraints of finite sky 
coverage and computation power. Instead, one divides the accessible £-range constrained by 
the sky coverage and the observing beam size into several bands {b}, and then estimates the 
band power Cf,, i.e., one approximates Ce in the form 

Q-aCf, (2-11) 
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where Cf^ is a chosen shape function characterizing the scale dependence in each band. For 
example, one can choose 

Cf = J^y (2-12) 

which leads to a scale-invariant form in each i band, i.e., i{i + 1)C^ = const V£ G b. With 
the approximation (2-11), one can rewrite equation (2-8) as 

Cspp' ~ / ^Cbf^pp'bjCi ,Bp,Bp'], (2-13) 

b 

where 



l^pp'b[Ci\ Bp, Bpi] — 2_^ 2_^ ^'< 



sh ^ 



i£b m=-i 

Bp£mB*,i^Yim{'^p)ye*mi'^p')- (2-14) 

If the beam is symmetric, then one has from equation (2-10) or (2-13) that 

Cspp' ~ 2_^ CbJCb[9ppr] CI \ -B^], (2-15) 

b 

where 

lCi,[9ppi] C| , -B^] = 
T.,eb'-^CtBjP,{cosepp,). (2-16) 

In the analysis procedure outlined above, the first problem arises in equations (2-2) and 
(2-3). Strictly speaking, what is convolved in reality is not the pixel temperature in Sp itself 
but the CMB signal in the time-stream 7^, i.e., 

00 I 

7i = X] X] BumaimYtmiy^t), (2-17) 

1=0 ■m=-e 

where Btim is the multipole expansion of the time-stream beam -Bt(x). This means that the 
experiment gives us a beam which moves on the sky as a function of time, t, and indeed may 
observe a different signal within the same pixel, p, depending on the orientation of the beam 
and the location of its center. We thus make a map which may have many different beams 
contributing to a single pixel. However, in our analysis formalism we must actually express 
this map as in (2-3), an observation of the sky with only a single pixel beam, Bp. Hence, for 
the Ci estimation, we need to find a way to estimate the pixel-beam expansion Bpim from 
the -Btfm, and this will be the focus of sections 4 and 8. 
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The second problem appears in equation (2-13). If the beam is not symmetric, the 
summation over m and the dependence on the pixel pair make the exact computation pro- 
hibitively expensive. To resolve this problem, in section 5 we introduce the pixel-pixel beam 
expansion B'^/g., which provides a consistent way to symmetrize asymmetric beams. This 
^i(cS) then replaces the Bj in equation (2-15), so as to approximate equation (2-13). 

On general grounds, the size of the observing beam is so small that, when necessary, 
we shall use the fiat sky approximation under the small- field limit. This means that when 
the size of a spherical patch is sufficiently small, the expansion of the beam in spherical 
harmonics is equivalent to a Fourier transform on a flat two-dimensional patch, i.e., 

5,™ = J dnB{^)Yem{^) = J dx^B{^)e-'''-^ = B{k), (2-18) 

and 

i = k=\k\. (2-19) 

Throughout the paper, we shall use a 'tilde' to denote the Fourier transform of a quantity. 



3. THE CRITERIA FOR BEAM SYMMETRY 

It is important to clearly deflne the level of asymmetry of an antenna beam. Consider 
the multipole expansion B^m of the beam. For a given i, the variance of B^jn about its mean 
over m is 

where Bl. is the mean of squares over m: 

^£(ms) = {\Bim\'^), (3-2) 

and -B£(gjni is the square of the mean over m: 

Bi(sm) = [Be{m)] = (Blm)'^. (3-3) 



Here -Bf(m) is the mean of B£m over m, and therefore can be either positive or negative. 
We also note that B"^,. is the power spectrum of the asymmetric beam, and that B'j,^ is 
equivalent to the power spectrum of a symmetric beam that is azimuthally averaged in the 
real space. Based on this, one can deflne an 'index of asymmetry' (lOA) as 






Be(ms) 



R2 

-°^(sm) 

£(ms) 



1/2 

(3-4) 



We see that wi varies from zero to one — the larger the w^, the more asymmetric the beam. 
We also note that if the beam is symmetric, then w^ is exactly zero. Thus for a given beam, 
W(^ provides us an objective measure of the level of its asymmetry. 

In certain situations, we need to combine data from two or more photometers with 
different beam shapes. We shall use a subscript i (i = 0, 1, 2, etc.) to denote the quantities 
obtained from different photometers. As an analog to equation (3-4), it proves useful to 
define an 'index of combined asymmetry' (lOCA) for all the beams as 



W, 



R2 

-°S^(sm) 

"R2 

-"S^(ms) 



1/2 



(3-5) 



where 



r2 

-°E£(sm) 



r2 



Z 



CiBie(tn) 



EC. 



r^2 

i£(ms) 



(3-6) 



(3-7) 



the Bif(m\ and B"^, 



i£(m) 



it.(vas) 



are the -B^(ni) and Bl. of photometer i respectively, 

ti(obs)/NET2 



E, [t,(obs)/NET^] ' 



(3- 



t 



is the total observation time of photometer i, and NETj is its noise equivalent tem- 



i(obs) 

perature (NET). Here B^^,s is the square of the noise-weighted mean of Bi, and B'^^,. 
is the noise-weighted mean of the squares of Bi assuming all the Bi are fully correlated. As 
one can see, the We varies between zero and one — the larger the We, the more asymmetric a 
Cj-weighted combined beam can be (depending on the detailed orientations of the beams in 
the temporal samples; we shall discuss this later). This also means that the lOA (ccj^) of an 
average beam with a weight Q for each Bi is always equal to or smaller than We, although 
the individual wn of Bi may be larger than We- If We = 0, then we know that all the beams 
are symmetric {wie = 0), and vice versa. 

For the purpose of power spectrum estimation, one can employ Wie (or We when combin- 
ing data of different observing beams) to decide if a simple symmetric-beam approximation 
is sufficient. For example, at £'s where zue ~ 0, we expect equation (2-15) to be adequate. 
On the other hand, at £'s where we (or We) deviates significantly from zero, one may need 
to employ equation (2-13). We shall further discuss these situations, and the use of the lOA 
(zue) and the lOCA (We) later. 
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4. THE AVERAGE PIXEL-BEAM EXPANSION 

4.1. The pixel-beam expansion 

We first estimate tlie 'pixel-beam expansion', Bp^^n-, from given observing beams, -B«m- 
A naive way to investigate this is to substitute equations (2-3) and (2-17) into the model 
(2-2), leading to 

AtpBp£mYem{^p) = BumXl.m\^t)- (4-1) 

This equation holds if and only if there exists one Xp for every x^ such that x^ = Xp. In 
this case, we have Bpim = Bum- This is of course true when the pixel size is infinitesimal, 
but is unlikely to be fulfilled in reality. Nevertheless, equation (4-1) is just the result of the 
modeling and therefore not necessarily a requirement in practice. In our formalism for the 
power spectrum estimation, the Sp is an unknown quantity to be estimated by using equation 
(2-5), so the actual relation between Bp^m and Bum should be also obtained through the same 
process. First, we substitute equation (2-1) into (2-5), and the CMB signal part yields 

Sp = CjvA^iV-Si, (4-2) 

where Cjq = C^ppi as defined in equation (2-9). Further substituting equations (2-3) and 
(2-17) into this result, we obtain 

BptmYtmi^p) = CnA N BumYtrni^t)- (4-3) 

This equation is completely general, and should be in principle satisfied when one tries to find 
the Bpf^rn from the given -Btfm- We thus see that equation (4-1) is just one of the solutions 
to equation (4-3), but not necessarily a requirement for the purpose of power spectrum 
estimation. 

In most cases, the noise rit in each temporal measure is nearly independent from the 
others, so the time-time noise correlation matrix Nu' is diagonal, with the tt elements equal 
to the noise variance at each time sample, i.e., 

Nu'=lil5{t-t'), (4-4) 

where fit is the standard deviation of time sample t, and 5{t — t') is a Dirac Delta. This 
allows us to simplify equation (4-3) as 

Bp£mYl.m\?^p) = / ^jtBumYimij^t)) (4-5) 

tep 

where ^t is the noise-estimated statistical weight at t: 

6 = ^^^. (4-6) 
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For simplicity, we shall take this white-noise assumption for further investigation. We con- 
sider the more general case of correlated noise in the Appendix, and show that this white- 
noise approximation is appropriate in most practical cases. The conditions for the use of 
this white- noise assumption will be also derived in the Appendix (see eq. [AlO]). 

To further simplify equation (4-5), we assume that xt = Xp Vt G p (i.e., the temporal 
measure 7t is thought of as a 'sample' of the pixel temperature Sp] see eq. [4-2]), so that the 
pixel-beam expansion can now be obtained as 

Bpim = 2_^ itBum- (4-7) 

The assumption, xj = Xp, for achieving this result will be relaxed in section 8, where we 
show that only an extra correction is required. 



4.2. The average pixel-beam expansion 

As will be shown, it proves useful to remove the pixel dependence of -Bp£m in the for- 
malism of the C^ estimation. We thus consider the noise-weighted average of -Bp£m over all 
pixels (c.f. eq. [2-5]): 

Bpim = H{U C^ Bpim), (4-8) 



where f/ = f/p is a contraction vector with entries all equal to unity, and 

H = {U^C],'U)-\ (4-9) 



We shall call Bp£m the 'average pixel-beam expansion'. We note that the subscript p in Bpim 
does not mean the pixel dependence as in the usual convention, but indicates that this is a 
mean taken over all pixels. 

With the white-noise assumption (eq. [4-4]), the Bpim can be calculated explicitly by 
substituting equation (4-7) into equation (4-8): 

Bpim = / ^XtBum, (4-10) 

t 

where 

Xt = ^^. (4-11) 

If the data are from a single photometer with a constant noise level, then equation (4-10) 
reduces to a simple linear average of all time-stream beams. If the data are combined from 
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different photometers, then the fit can be approximated as (c.f. eq. [3-8]) 

where NET^ is the NET of the corresponding photometer at time t, and 5t(obs) is the inte- 
gration time of the temporal observation at t. If the integration time remains unchanged 
among photometers, then the Ht in equation (4-11) can be simply taken as the NET of the 
corresponding photometer. We also note that with the definition (4-12), equations (3-8) and 
(4-11) can be related as 

= ^Xt, (4-13) 

tei 

meaning that Q is the total noise-estimated weight of photometer i. 

We note that in cases where both the shape of the experimental beam and its orientation 
relative to the pixel are roughly constant throughout the observation, we have a reasonable 
approximation (see eq. [4-10]): 

Bp£m ~ Bum- (4-14) 

In other cases, equation (4-10) will need to be employed, for example, when the relative 
orientation between the asymmetric beam and the pixels changes, or when data from different 
photometers are combined together. We also note that even if all the beams Bi of different 
photometers are symmetric (i.e., Wie = Wi = 0), the Bpfm may still have pixel dependence 
due to the various relative contribution of Bi within different pixels (see eq. [4-7]). In such 
cases, one will need to consider equation (4-10), and a simple formalism like equation (2- 
15) will be invalid for the estimation of the CMB angular power spectrum, since the 5L is 
different on each pixel. As will be shown, the formalism we shall develop is also capable of 
dealing with this situation. 



4.3. Useful Limits 

We now derive useful constraints on the magnitude of the average pixel-beam expansion 
Bpim- In the small-field limit, the power spectrum of Bpim can be written as (see eqs. [2-18], 
[3-2], and [4-10]) 



Bpi^nis) = - I dip 



Tl 



Y.^tBt{\^) 



(4-15) 



where ip is the phase angle of k on the ring |k| = k. We first consider single-photometer 
experiments. In this case, if the beam pattern remains the same throughout the entire 



-12- 



observation but with only different orientations at different t, then we can rewrite Bt as 

Bt = A{l3t)Bo, (4-16) 

where A{f3t) is the rotation matrix, (3t is the rotation angle at time t with respect to t = 0, 
and Bq is the shape of the time-stream beam at t = 0. Substituting this into equation (4-15) 
gives 



B 



p£{nis) 



- I dip 

TV 



27r 



f{P)A{P)Bo{k)dp 



(4-17) 



p27r 



where /(/9) is the weighting function of a rotation angle P, and satisfies Jq"" dPf{P) = 1. It 
is then straightforward to show that the function f{f3) that minimizes the right hand side of 
the above equation is f{P) = 1/2-11, leading to 

1 2 



p£(ms) 



dipBoik) 



r2 

-D0^(sm)5 



(4-18) 



where Bq^,^s is as defined in equation (3-3). On the other hand, the function /(/5) that 
maximizes the right hand side of equation (4-17) is /(/9) = 5{I3 — f3o) (Dirac Delta, f3o G 
{0, 2n}), and this gives 



B 



pl(ins) 



n 



dip 







5o(k) 



R2 

-°0^(ms)5 



(4-19) 



where B^^^,-. is as defined in equation (3-2). These results tell us that when the pixels are 

2 

scanned almost uniformly in all directions, then the resulting -Bp^fms) should be closer to 



B 



p£(rns) 



R2 

-°0^(sm) 



p£(ms) 

. When the pixels are scanned with an almost fixed direction, then the 



resulting B^^,^^-. should be closer to B^^,^^-, 



B. 



o£(ms)- Thus, we have a good check of 



the numerically calculated Bp£m from equation (4-10) (or eq. [4-17]), i.e., a constraint on the 

;2 

p^(rns) ' 



2 

amplitude of B 



B, 



0£(ms) 



> B 



p^(ms) 



>B' 



(sm) ' 



or equivalently. 



1 > 



B 



p£(ms) 

"^2 

-^O^(ms) 



> 1 



ZUf, 



(4-20) 



(4-21) 



where w^i is the lOA of Bq. For symmetric beams, all the equality signs hold. In experiments, 
one can take Bq as the measured beam, and then use equation (3-4) to calculate Wq£. 

When we combine data from two or more photometers with different beam shapes, 
following the same line of development as above gives (see eqs. [3-6], [3-7], [4-10], and [4-13]) 

:r2 



-^S£(ms) — Bp£{ins) — -DS^(sm)5 



(4-22) 



or equivalently, 
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-r2 

1 > -^f^ > 1 - Wj, 



(4-23) 



where We is the lOCA defined in equation (3-5). We shall further discuss the use of these 
limits later. 



THE PIXEL-PIXEL BEAM EXPANSION 



5.1. Formalism 

In the data analysis procedure briefiy demonstrated in section 2, the effect of asymmetric 
beam convolution manifests itself in equation (2-8). However, the summation over m and 
the dependence on the pixel pair make it computationally expensive. Therefore, we prefer 
to use the form of equation (2-10) as an approximation. This can be achieved by replacing 
the Bj in equation (2-10) with a 'pixel-pixel beam expansion' B'j,^-., which we shall derive 
in this section. 



First, one can replace the Bf in equation (2-10) with 

_ '^'^YjrnBpimBl,^^Yi>m{^p)Y^^X'^p, 



B^ 



pp'i 



{2i+ l)Pi{cose, 



(5-1) 



pp 



so that equation (2-10) is equivalent to equation (2-8). In the small-field limit, equation 
(5-1) becomes 



J 



B 



pp'i 



B, 



kAx,<fo;Bl„,{k) 



pp'k 



Jo{kAx) 



(5-2) 



where 



J 



kAx,ipo;Bl,{k) 



^^dip\?R: i32p,(k) cos[A;Aa;cos(v?-v?o)] 



-3 



i3p ,(k) sin[A;Aa;cos(v9 - v^o)] 



(5-3) 



Ax = Xp — Xp/, Ax = |Ax| = 6ppr, ipo is the phase angle of Ax, Jq is the Bessel function of 
the first kind of integral order 0, Bpp,(k) = Bp(k)B*,(k), and 3fJ and Q' indicate the real and 
imaginary parts of B^p, respectively. We notice that J'[kAx,ipo; 1] = Jo{kAx). Therefore if 
the beam is circularly symmetric and remains the same on all pixels, i.e., B^ ,{\<.) = i?|, then 
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J[kAx, (po; i3pp,(k)] = J[kAx, ipo; Bl] = Jo{kAx)Bl, so that -8^^/;. in equation (5-2) becomes 
Bl exactly as required. 

To save computation time and memory when estimating Ci, we need to remove the 
dependence of -B^p/^ on the particular choice of a pixel pair (xp,Xp/). We achieve this by 
taking the average of B^^,^ over all possible (xp, Xp/) pairs: 

Be{eS) = {Bpp'e) ■ (5-4) 

We call this -B^/^ff) ^^^ 'pixel-pixel beam expansion'. 

Even with this, equation (5-4) together with equation (5-2) is still computationally 
expensive and may not be feasible. Therefore we further simplify the formalism in the 
following way. First, we remove the dependence of Bpp,(k) in equation (5-2) on pixel pairs, 
by replacing it with a noise- weighted average (c.f. eqs. [2-5] and [4-8]) 



^l' = {U^C^'UU^C^'U)-\U^C^'Bl,C^'U). (5-5) 



Here the subscript pp' in 8^^, does not mean the pixel pair dependence as in the usual 
convention, but indicates that the mean is taken over all pixel pairs. With this replacement, 
equation (5-2) is now only a function of Ax for a given i = k. Thus when evaluating equation 
(5-4), we can classify all possible Ax into several groups of different Ax, each with several 
subgroups of different ipo. This gives 



?2 — \:^ 



J 



kAx,^o;l3l,ik) 



^e{cs) = 2^ 9{^x,ipo) 7777X3 ' (^-6) 

Ax,ipo 



JoikAx 



where g{Ax,ipQ) is the weight of the configuration (Aa;,(/9o), i.e., the number of pixel pairs 
with Ax and v^o; divided by the total number of pixel pairs. It satisfies X^Ax w-, di^^^ Vo) = 1- 
This algorithm can normally reduce the number of operations in equation (5-4) by several 
orders of magnitude, because the element number of {(Ax,(^o)} is normally several orders 
below that of {(xp, Xp/)}. In addition, if the number of pixels is large enough as in most cases, 
then ipo is nearly uniformly distributed between and 27r for every given Ax, depending on 
the relative locations of all pixels. In this case, after the summation over ipQ at each given Ax 
in equation (5-6), the first term inside the integral in equation (5-3) (which enters eq. [5-6]) 
becomes 3f? i3p ,(k) Jo (A; Ax) and the second term vanishes. Thus the Bessel function in 
equation (5-6) can be removed and we have 

1 



i?.V) ^ ^ y„ ^ 



^L(k) d^. (5-7) 
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With careful simplification of the real part of equation (5-5), we also find that 

7 2 



3? 



BL,(k) 



pp 



Bp{k) , (5- 



where -Bp(k) = Bpim as defined in equation (4-8). We note that the average over all pixel 
pairs (the left-hand side of eq. [5-8]) is now reduced to the average over all pixels (the 
right-hand side). This further enables us to simplify equation (5-7) as 



r2 _ 1 



TX 



5p(k) dif = 5p^(,n,), (5-9) 



— 2 

where the last step uses the definition (3-2), and the -Bp£(ms) is readily evaluated in equation 

— 2 

(4-15). When calculating Bp^/^^^-., one can take the form of equation (4-17) to save compu- 
tation time. We note that the approximation sign above will become equality when ipQ is 
uniformly distributed between and 27r. In section 9, we shall numerically verify this result. 

With such, now we can use the form of equation (2-10) to approximate equation (2-8) in 
the presence of asymmetric beams or when combining data with different symmetric beams. 
In other words, we have equation (2-8) being approximated as 

Cspp' ^ 22 ~A CiB\,Q)P^{cOS Opjy). (5-10) 



1=0 



Att 



Furthermore, as illustrated in equations (2-11) through (2-16) and the context, one normally 
divides the i range under investigation into several bands, due to the finite sizes of the sky 
coverage and the observing beam, as well as the limited computation power. Using this 
formalism, we can approximate equation (2-13) using equation (2-15) with its -B| replaced 
by the B^.^^^. calculated above. This gives 

Cspp' ~ 2_^ CblCb[6pp'; C| , -B^(cff)]- (5-11) 



5.2. Uncertainties 

When making the approximation (5-10), we inevitably induce errors in the basis B'j,Q^Pi{cos Oppi 
for each pixel pair. These errors can be represented as 

^ = 1 ± a,, (5-12) 



-^^(eflf) 
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where ai is the normalized standard deviation of the errors. This deviation can be simulta- 
neously evaluated while one performs equation (5-6), i.e., 



al 



J 



kAx, ipo] B: 



Since Ci appears in combination with B'^,g-.Pi{cos9pp') (see eq. [5-10]), we know that ai 
basically quantifies the bias in Ce for each individual pixel pair. Nevertheless, the resulting 
bias in the final Ce estimates by using the approximation (5-10) together with the likelihood 
analysis (see eq. [2-6] and context) may be much smaller than a^, because the resulting C^ 
is a consequence of the contribution from all pixel pairs. For example, if all pixel pairs 
contribute to the likelihood function (2-6) as a linear combination of B^ ,fPi{cos6ppr), then 
the resulting bias in Ce will be as small as ae/J^p, where A/^ is the total number of pixels. 
Although we know that reality is not like such a simple case, we can still quantify the bias 
of approximation (5-10) using numerical simulations. 

Similarly, we can consider the errors in the band power Cb for each individual pixel pair, 
resulting from the approximation (5-11). Since Cb is coupled with ICb (eq. [5-11]) or /Cpp'f, 
(eq. [2-13]), the errors in Cb for each individual pixel pair may be quantified by comparing 
/Cfe and ICppfb, as we did for B'^.g. and -B^p/^. However, as argued earlier, the result calculated 
in this way quantifies only the errors in Cb for each individual pixel pair, and the real bias 
of the approximation (5-11) together with the likelihood analysis may be much smaller. We 
shall quantify the real systematic bias of this approximation in section 9, using numerical 
simulations. 



6. SYMMETRY VS. ASYMMETRY 

In this section, we investigate the conditions under which one needs to employ the 
formalism for treating asymmetric beams, i.e. the formalism we developed in the previous 
two sections. We first consider the case where the data to be analyzed is from only one 
photometer. From equation (5-9), we know that the pixel-pixel beam expansion -B^^^ff) '^^^ 

2 

be approximated by -Bp£(ins)? ^^^ therefore should be also constrained by equation (4-21), 
resulting in 

l>#^>l--oV (6-1) 



-^0£(ms) 



This implies that if we simply use -Bo£(ms) (where Bq is the measured beam shape from the 
experiment) as the -B|(eff) ^^ '-'^^ formalism, then B'j,^^-. will be overestimated by at most 
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Furthermore, we consider the errors in Ci resulting from this effect. In our formahsm, 
the beam convolution appears as the multiplication of B'jr^s^ and C^ (see eq. [5-10]), so the 
errors in C^ can be expressed as 

Taking dBz'^g-. = Bq3. — B^/^^-., we have from equations (6-1) and (6-2) that 

0>6c,> -wl,. (6-3) 

This means that when we use B?.f, . as the B], g. in our formalism, then the resulting Ce 



O^(ms) ^° "^^^ -^e{eS) 
2 

eflf) 

2 _ l-wl 

O^(mid) — 7^32 I R-2 ~ 1 _ „2 /9-^0^(ms); 



at a given I will be underestimated by at most tUg^. To share this error on both sides of a 
mis-estimated C^, we can choose the Bj, rr-. to be 



r2 _ f _ ^ '^0£ p2 /(- /,N 

-^O^fmid) — p-2 I R-2 ~ 1 _ „2 /9-°0<?(ms)) l^'^J 



SO that the resulting error in the C^ estimates is now constrained as 

2 2 



'W 



2 - 2n7Q^ 



> <5c,(mid) > -^. (6-5) 



When Wq^ <^ 1, we have |(5c^(mid)| ^ '^m/^- If the beam is symmetric, then all the equality 
signs above hold and Bl^^^^^ = B^^^^,^^ = B^j^^^^^ = B\^^y 

Following the same line of logic, we now consider the cases where the data to be analyzed 
is combined from two or more photometers. In this case, it is also straightforward to show 
that if we choose the B'j,^-. to be (see also eqs. [3-6] and [3-7] for definitions) 

^E.(.nid) - ^2 ^2 - 1 _ -^2/2^S^(ms), (^-6) 

L£(ms) L£(sm) t ' 

then the errors in the C^ estimates are constrained as 

> ^^(mid) > TT- (6-7) 



2 - 2W| ^'^"'''^^ - 2 

When >V| <^ 1, we have |(5c^(mid)| ^ ^J/2- If all the beams are symmetric {wu = 0), then 
all the equality signs above hold. 

As a result, we see that if the cc7q^/2 (or VV^/2) is well below the tolerated maximum error 
of Ci, then we can use -Bo£(mid) ('^^ -^l^(mid)) ^^ ^^^ ^e(cs) ^^ ^^^ symmetric-beam formalism, 
i.e., we can simply use equation (2-15) with Bj = -Bo£(mid) (°^ -^l^fmid))' without the need 
of going through the procedure developed in sections 4 and 5. The associated errors in the 
final Ci estimates will be constrained by equation (6-5) (or eq. [6-7]). 
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7. UNCERTAINTIES FROM BEAM MEASUREMENT 

It is inevitable for any experiment that there are uncertainties in the measurement 
of the beam shapes. It is therefore crucial to quantify the uncertainties in the final Ci 
estimates resulting from this beam shape uncertainties. For a given beam -B(x), consider an 
uncertainty e in the full width at half maxima (FWHM), and assume that the uncertainties 
at all iso- height contours of the beam are a fixed fraction of the contour sizes, i.e., 

This uncertainty in the beam shape will then be transfered to the multipole space as the 
uncertainty in £ at a given height Bfm'- 



This results in the uncertainty in Bf at a given d, 

^B^ = -5T = —52 1- (7-3) 



We then consider the change in the C^ estimates: 

dC 



A d(^i ,„ ,^ 

Ac, = ^. (7-4) 



Since the beam convolution occurs as the multiplication of -B| and Cn (see eqs. [2-8] and 
[5-10]; here we have dropped the subscript '(eff)' for concise notation), we know that the 
resulting uncertainty in C^ is 



^Q 



R2 
52 



- 1- (7-5) 



This means that if the beam size is mis-estimated by e (i.e., the actual size is 1 + e times the 
measured size), then the resulting C^ estimates will be 1 + Aq times the real C^. Thus for 
a given uncertainty in the beam measurement e, one can employ equation (7-5) to estimate 
the resulting uncertainty in the final C(, estimates. We also note that the banding of I does 
not affect this result, as we shall show in section 9.3. 

We note from the above result that the type of uncertainty we assumed in the measure- 
ment of the beam size (eq. [7-1]) induces an error that is correlated between the Ci estimates 
in all 1 bins, although the magnitude of the error depends on l. If the power of the beam 
B\ monotonically decreases with £, a measured beam size slightly larger than the real value 
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will produce larger Ci estimates at all I bins. Although this is not the most general kind of 
beam size error, it is quite common. 

We now investigate certain special cases. In situations where dBj/dd is not changing 
much within M, i.e., 

M « !^0±l)£ (7-6) 

we can approximate equation (7-3) and thus (7-5) as 

where equation (7-2) has been employed. For a symmetric Gaussian beam -B(x) = exp{—x'^/2g'^) 
equation (7-7) becomes 

Ac.(G) ^ ^h(G) = -2ef^'^^ (7-8) 

while the condition (7-6), for |e| ^ 1, leads straightforwardly to 

1 



< iU = ^ . (7-9) 



Here we have again used the small-field limit. When combined with the condition (7-9), 
we find that approximation (7-8) breaks down when |A^ .^J is comparable with unity. In 
particular, we investigate the accuracy of approximation (7-8), by comparing it with equation 
(7-5). We find for |e| < 20% that the approximation is accurate within 10% error if 

^<^(G:io%) = (0.44 + 0.8|6|)£^G), (7-10) 

where a,*,^-. is given in equation (7-9). For example, if e = 10% and the Gaussian beam has a 
FWHM of 10 arcminutes (i.e., q = 1.24 x 10~^ radians), then approximation (7-8) is accurate 
within 10% error when i < ^*q.iq%) ~ 940. Under the condition (7-10), one can see from 
equation (7-8) that, for an approximately Gaussian beam, the resulting uncertainty in the 
final Ce estimates increases in proportion to the uncertainty in the beam measurement e, the 
square of the beam size g^, and the square of the multipole number i"^. 



8. DECONVOLUTION OF THE PIXEL SMOOTHING 

We have not dealt with the smoothing effects due to the pixelization of the map, when 
translating the data from the temporal to the pixel domain (see eqs. [2-5] and [4-3]). Because 
convolving a CMB map with a Dirac Delta (5(x — xi) will shift the original temperature at 
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X to a new location x + xi, we know that Yfm{'^) = 5^m(xi)Y£m(x + xi) where Sem{^i) is the 
multipole expansion of 5(x — xi). This allows us to rewrite equation (4-3) as 

Bp£m = CnA N Bum^lm{^p " X(). (8-1) 

Substituting this into equation (4-8), we obtain 

Bp,m = HU^C],'diB.g{CNA^N-'Bs,J, (8-2) 

where Bsem = Btem^imi^p — ^t) is a Aft by Afp matrix, and diag(M) is a vector whose entries 
are the diagonal elements of the matrix M. These results are completely general. Without 
further information about A^^^ or Bsim, equation (8-2) can not be simplified, mainly due to 
the involvement of 6em{'^p — xt). 

With the white-noise assumption (see sec. 4.1), we have equation (8-1) simplified as 



t&p 



and equation (8-2) as 



Bupim — 2_^ XtBurn^tmi^p^t — Xt), (8-4) 

t 

where x^g^ is the central coordinates of the pixel p that covers x^, and ^t and Xt are as defined 

in equations (4-6) and (4-11) respectively. Here we use the subscript '11' to distinguish these 

results from those in equations (4-7) and (4-10). In the real space, equation (8-4) is equivalent 

to 

^np(x) = ^ Xi5i(x - Xp3t + Xi), (8-5) 

t 

meaning that Bnp is the noise-weighted average over the time-stream beams Bt that are 
shifted by Xpgj — x^ at each time t. This implies that our formalism developed previously 
is still available, requiring only a modification that takes into account the detailed locations 
of the temporal hits x^ with respect to the pixel centers Xpg^. Thus we have relaxed the 
assumption x^ = Xpg^ that was made to achieve equations (4-7) and (4-10). 

In most cases, both Mp and Mt are large, and the beam shape Bt of each photometer does 
not change much within several successive pixels. This results in the fact that in determining 
the Bnp in equation (8-5), each beam configuration A(/9)-Bo(x) (see eq. [4-16]) appears at 
a set of Xt which have offsets Xpg^ — x^ distributed within a region that is confined by the 
pixel shapes. If all pixels have the same shape, then this is equivalent to convolving each 
A{I3)Bq{:x) with a top-hat like window whose boundary is defined by the pixel shape. As a 
result, we can approximate equation (8-4) as 

Buplm ~ BpijyJIirn, (8-6) 
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where Bp^m is as defined in equation (4-10), and Il^m is the multipole expansion of 

n(x) = 5^5(xp3i-Xi). (8-7) 

t 

Tlie same also applies to the simple case where all time-stream beams B^ are the same. We 
thus see that the Hem in equation (8-6) serves as an extra convolution (apart from the time- 
stream beam convolution) of the CMB signal due to the pixelization of the map. With such, 
we can now easily incorporate this extra smoothing effects into our formalism by replacing 
our B^^f^. with (see also eq. [5-9]) 

-^n^(eflf) ~ -Snp£(ms) (8"8) 



^2 



-^p^(ms)n^(ms) (8-9) 

Be{eS)^e{ms)y (8-10) 



-ttS 



where -Bnp£(ms) is the power spectrum of the Bjip defined in equation (8-5), and Tiji^g\ is the 
power spectrum of the n(x) defined in equation (8-7). We note that in the limiting case 
where x^ = Xpgj, we have Hpim equal to unity for all i and m (since the multipole transform 
of (5(x) is unity), so the smoothing effects disappear, and we have exactly -Bn£(efj) = -B^(cff) 
(see eq. [8-10]). If the pixels do not have exactly the same shape, as in the case on any large 
patch of the sphere (e.g., pixelized by HEALPix, Gorski et al. 1999, or by Igloo, Crittenden & 
Turok 1998), then we can use equation (8-8) together with equation (8-5) to obtain -Bn^(cff)- 
If the pixel beam or the pixel shape remains roughly the same for all pixels, then we can use 
equation (8-9) together with equations (5-9) and (8-7) to calculate -Bn£(cff)- 

If all the pixels have the same shape which is a regular square of size q in radians, then 

we have 

2 .. 8 ['- ^^ sin^ [(£^cos0) /2] sin^ [(£^sin0) /2] 

^ ' 71 Jo i^q^lcos^ (p)[sm (p) 

An accurate approximation to this result is 

(£^2.04 



n£(ms)(^) ~exp 



18.1 



[1 - 2.72 X 10-2(£^)2] . (8-12) 



The accuracy of this fit is within 0.3% error for i < lAn/q. For example, if the pixel size is 
5x5 square arcminutes (i.e., ^ = 5 arcminutes ~ 1.45 x 10^^ radians), then the above fit is 
at 99.7% accuracy for £ < 3024. 
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9. NUMERICAL VERIFICATIONS 

9.1. The pixel-beam expansion 

In this and the following three subsections, we will employ an elliptic Gaussian beam 
with a short-axis FWHM of 5 arc-minutes and a long-axis FWHM of 20 arc-minutes, to 
demonstrate certain key points developed previously. We first investigate the pixel-beam 
expansion of a given pixel resulting from different scanning strategies, i.e., to investigate the 
dependence of the pixel-beam power spectrum (4-17) on the function f{f3), and to verify 
the results given in equation (4-20). We note that although those results are given for 
the average pixel-beam expansion, we expect the pixel-beam expansion to carry the same 
property since equation (4-7) has exactly the same form as equation (4-10). Figure 1 shows 
two different configurations of beam scanning on a given pixel. In case A, the pixel was 
hit twice by the same beam pattern, but with different orientations of a separation angle 
a. That is /(/3) = [S{P) + 6{P — a)]/2. We shall investigate the cases a = 15, 45, and 90 
degrees. In case B, the pixel was hit evenly in four different directions. That is /(/3) = 
[6{f3) + 6{f3 - 45°) + 6{f3 - 90°) + 6{f3 - 135°)]/4. 

Figure 2 shows the lOA of the pixel beams in cases A and B, as defined in equation 
(3-4). As one can see, the pixel beam has the largest asymmetry (largest w^) when the pixel 
is hit by a beam with only one direction (the dashed line). When the pixel is hit by beams 
of two different directions (case A in Figure 1 ), the asymmetry decreases {we decreases) if 
the separation angle of the two directions a is closer to 90 degrees (see the dotted lines in 
Figure2 ). When the pixel is scanned with four different directions (case B in Figure 1 ), the 
resulting effective beam is nearly symmetric {w£ ~ 0) up to £ ~ 1000, and has the lowest 
level of asymmetry (the smallest w^) . 

Figure 3 shows the power spectra of the pixel-beam expansions with different scanning 
strategies. As we can see, the power spectrum of the pixel-beam expansion has a maxi- 
mum given by equation (3-2) (see also eq. [4-19]), when the pixel was scanned with only 
one direction. On the other hand, the power spectrum of the pixel-beam expansion has a 
minimum given by equation (3-3) (see also eq. [4-18]), when the pixel was scanned evenly in 
all directions (note that the dot-dashed line in Figure 3 almost coincides with the solid line). 
This verifies our results given in equation (4-20). By comparing Figure 3 with Figure 2, 
we also learn that there is a strong correlation between we and the -B| of a pixel — when 
the pixel is scanned by a same beam pattern with more different directions, the level of the 
effective beam asymmetry (we) decreases, and so does the power spectrum of the pixel-beam 
expansion (Bj). 

We also note that according to equation (6-5), the lOA of the original time-stream 
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Fig. 1. — Beam configurations on a given pixel, resulting from different scanning strategies. 
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Fig. 2. — Indices of asymmetry of the pixel-beam expansions, as functions of the multipole 
number L 
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beam Bo (the dashed hne in Fig. 2) also tells us beyond what i we need to worry about 
the asymmetry of the beam. For example, when i < 600, we see that wi < 0.3, giving 
ro|/2 < 0.05. This means that if we simply use the B^^^^.^^ = 2/[5q"^^^^) +5g-|^)] (eq. [6-4]) 
as the B'^,g. in the formalism (5-10), then the maxima error in the final C^ estimates is 
guaranteed to be within about ±5% for i < 600. 



9.2. The pixel-pixel beam expansion 

We now use the elliptic Gaussian beam of 5 by 20 arcminutes to verify some important 
results in section 5 — mainly equation (5-9). Consider a square map of size 10° x 10°, with 
a square pixel size of 5 arcminutes. Referring to equation (5-6) with such a map. Figure 4 
shows how (fo is distributed at each Ax. In the figure, each dot labels the {Ax, ipo) that is 
sampled by the map. As one can see, (fo is nearly uniformly distributed for any given Ax, 
except when Ax is close to the boundaries constrained by the pixel and field sizes. Because 
of this nearly uniform distribution, we achieved equation (5-9) from equation (5-6). 

More precisely, we carried out equation (5-6) to obtain Bl^^s, and calculated the right 
hand side of equation (5-9) to obtain B'j,.. Here we have used the elliptic Gaussian beam 

directly as the -Bp(k) = B^^rn in the due calculations. We found that Bis agrees with 
^HcS) with more than 99% accuracy for i = 0-2000. 

We have also calculated the average deviation o"£ of B^ ,^ from Bl^^. for each individual 
pixel pair, using equation (5-13). The result is shown in Figure 5. First, we see many spikes 
in a£. This is due to the zeros of the Bessel function Jq, which appears at the bottom of 
equation (5-13). These spikes should be neglected, as in reality no such singularities appear 
in our analysis pipeline. We note that these spikes have the same origin as those presented 
in Hanany et al. (1998), where a similar situation was considered. Second, as addressed 
previously, although the ae obtained from equation (5-13) can be as large as comparable 
to unity, the real errors in the final Ce estimates by using the formalism (5-11) with the 
approximation (5-9) will be much smaller than this value. This is because the a^ here tells 
only the mean discrepancy of B'j,^-. for each individual pixel pair, and may average out when 
all pixel pairs come into account in the likelihood analysis. In section 9.5, we will numerically 
justify this and thus the accuracy of employing equation (5-11) with (5-9). 
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Fig. 3. — Power spectra -B| of the pixel-beam expansions with different scanning strategies. 
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Fig. 4. — Distribution of ipo as a function of Ax. Each dot labels the {Ax, ipo) that is sampled 
by a square map of 10° x 10°, with a pixel size of 5 arcminutes. 
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9.3. Uncertainties from beam measurement 

In this section, we will numerically verify the results in section 7. First, we use a 
symmetric Gaussian beam with a FWHM of 10 arcminutes, to investigate the approximation 
A^ ,^s given by (7-8), as a comparison to the exact result Aq(g) given by (7-5). Here we 
take the uncertainty in the beam measurement to be e = 10% (eq. [7-1]). As one can see 
in Figure 6, the approximation breaks down towards the limit i*^^-. given by equation (7-9). 
For i <^ £^Q^, the approximation A^ ,^-. reproduces the correct result Aq(g). By comparing 
A^ /Q^ and Aq(g), we calculate the 10% accuracy hmit ^(g.io%) (^^e dot-dashed line), at 
which A^ ,^JAci,(G) ^ 1 = 10%. In addition, by varying the value of e between ±20%, we 
obtain the result presented in equation (7-10). That is, for a symmetric Gaussian beam with 
an uncertainty of e in size, the approximation (7-8) for the resulting uncertainty in Ce is 
accurate within 10% error for i < ^(g.io%) = (0-44 + 0.8\e\)i*^y 

Now we investigate the case where the beam is asymmetric. We use an elliptic Gaussian 
beam, whose long- and short-axis FWHM's are 20 and 5 arcminutes respectively. This beam 
is first convolved onto a simulated CMB map of size 10° x 10°, with a pixel size of 10 
arcminutes. The underlying cosmology is an inflationary model with [Q]^, Qcdm,^ A, n,h) = 
(0.07, 0.61, 0.23, 1, 0.60), normalized to the CORE DMR. A random Gaussian noise of lOO/xK 
is then added into each pixel. We call this simulation (1). We repeat the same procedure 
again except that this time the beam size is increased by 10%, i.e., e = 10%, to obtain a 
simulation (2), where the CMB and noise realizations are exactly the same as those used in 
simulation (1). We then analyze both simulations using the procedure outlined in section 2, 
with the approximation (5-11). The resulting uncertainty in Cg can thus be calculated using 
equation (7-4) as 

Ac,,., = p^-l. (9-1) 

where the subscripts (1) and (2) indicate results from the two simulations. The results are 
shown as crosses in Figure 7. Also plotted is the result using equation (7-5) (the solid line), 
which we label with a subscript (b). It is obtained directly by varying the beam shape 
with e = 10%. As one can see, the crosses are highly consistent with the solid line. This 
means, flrst, that the asymmetry of the beam does not affect our result given by equation 
(7-5). Second, the banding of i does not affect the result, so we can use equation (7-5) 
as an estimate for the uncertainty in the band power Cb resulting from that in the beam 
measurement. This is also an important support to the fact that the banding of i does not 
affect the general relation Ci oc -B^^ (see eqs. [5-10] and [5-11]). We have also verified that 
the sizes of the error bars in the Ce estimates between simulations (1) and (2) do not change 
by more than 4% for i < 1200. Thus we know that when Aq is small, the uncertainty in the 
beam shape measurement does not affect the sizes of error bars significantly, but does affect 
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Fig. 5. — Mean discrepancy a^ of -B^p/^ from B'^^^g. for each individual pixel pair 
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Fig. 6. — Uncertainty in the Ce (solid line; given by eq. [7-5]), resulting from an uncertainty 
of e = 10% in the beam shape measurement of a Gaussian beam with a FWHM of 10 
arcminutes. Also plotted is the approximation (7-8) (dashed line). The dotted vertical line 
indicates the limit of the approximation £^qn given by equation (7-9), while the dot-dashed 
vertical line shows the 10% accuracy hmit ^(gio%)' which is well fitted by equation (7-10). 
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the amplitudes of the C^ estimates. On the other hand, when Aq is large (comparable to 
one), the signal to noise ratio may be affected and so may the error bar sizes. 



9.4. Deconvolution of the pixel smoothing 

We now test the formalism of deconvolving the smoothing effect due to the pixelization 
of a map. This is to verify equation (8-9), with equation (8-12) as an approximation in cases 
where the pixels are regular squares. We consider a square CMB map of size 10° x 10°, with 
regular-square pixels of size 10 arcminutes. We first simulate a time-stream of the CMB 
signal 7(, that is convolved with an elliptic Gaussian beam of 5 by 20 arcminutes in FWHM 
(same as the one used in previous sections). For each temporal sample, we then add Gaussian 
random white noise rit with 5% in RMS amplitude. In this run, we require the temporal 
samples to be exactly at the centers of each pixels, i.e., Mt = Mp and x^ = Xpg^, such that 
m = d (see eq. [2-5]). We call this simulation (0). In a second run, the procedure is the 
same except that the CMB temporal samples now have offsets with respect to the centers 
of each pixels, i.e.. A/* = A/^ but x^ ^ ^p3t with x^ — Xpg^ randomly distributed within a 
square of size 10 arcminutes. We call this simulation (1). In third, fourth, and fifth runs, the 
procedures are the same as simulation (1), except that the numbers of temporal samples in 
each pixels are now 3, 10, and 200 (i.e., Nt = SJ\fp, lOAp, and 200A/J,) respectively, instead of 
one. We denote these as simulations (3), (10), and (200) respectively. All these runs are then 
analyzed in the same way, using the procedure outlined in section 2, with the approximation 

9 

(5-11) and -B|(gm = -Bp^(jns) (eq- [5-9]). Therefore the ratio 



Ce{j) 



a 



m 



3 = 1, 3, 10, 200, (9-2) 



will quantify the smoothing effect due to the pixelization of the map. We plot this ratio in 
Figure 8, as a comparison to the n^/n^g) given in equation (8-12). 

As one can see and expect, the smoothing effect approaches the top-hat-window approx- 
imation when the number of temporal samples per pixel increases. When it is larger than 
10, as in most real situations, the top-hat-window approximation appears to be a good one. 
Also plotted at the bottom-left corner is the n(x) given by equation (8-7) for j = 10. The 
nearly uniform distribution of Xpgj — x^ shows the appropriateness of the top-hat-window 
approximation. Thus we have verified that the approximation (8-9), with equation (8-12) 
for cases where pixels are regular squares, is indeed a good approximation. For an obvious 
mathematical reason (see sec. 8), we know that this formalism can be further extrapolated 
for cases where pixels do not have regular shapes but the time-stream beams have roughly 
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Fig. 7. — Uncertainty in the C^ resulting from that in the beam measurement. The horizontal 
axis is the multipole number i. The crosses, AQ(a), are results based on simulations using 
equation (7-4) (or eq. [9-1]), while the sohd hue, Aci(h), is obtained directly from the beam 
shape using equation (7-5). An elliptic Gaussian beam, with long- and short-axis FWHM's 
of 20 and 5 arcminutes respectively, is used. The uncertainty in the beam shape is e = 10%. 
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Fig. 8. — Smoothing effect due to the pixelization of a CMB map. The ratios C£(j)/C^(o) are 
compared to the approximated smoothing window n^^j^g-,, with j = 1, 3, 10, 200 representing 
the number of temporal samples per pixel in different runs. The horizontal axis is the 
multipole number i. The square at the bottom-left corner shows the n(x) for j = 10. 



-30- 

the same shape. In such cases, equation (8-7) can be employed to obtain the ^^/^^g^ for the 
use of equation (8-9). 



9.5. The MAXIMA experiment 

In this section we demonstrate the application of our formalism using the data from the 
MAXIMA-1 experiment (HOO). Figure 9 shows the antenna patterns Bi {i = 1, 2, 3, 4) for the 
four photometers used in the analysis of the MAXIMA-1 data. Details of the measurements 
of these beam shapes are given in HOO. As one can see, the beams are more symmetric 
towards their centers. 

Figure 10 shows the level of asymmetry of the beams. The dotted lines are the lOA zun 
of each individual beam B^ (eq. [3-4]), and the solid line is the lOA of the noise- weighted 
combination of all of them, i.e., the we of the average pixel-beam expansion Bp£m (see eqs. [4- 
10], [4-11], and [4-12]). The dashed line is the lOCA Wi from all Bi (eq. [3-5]). Here the 
relative weight of each beam is Ci ^ C2 ^ Cs • C4 = 64 : 64 : 81 : 36 (see eq. [3-8] and HOO). 
As we can see, the beams are nearly symmetric (W^, zo£ ^ 0) at low i, but less so at larger 
i. At £ < 800, which is the range of i discussed in HOO, the asymmetry is less than 15%. 
The figure also confirms the fact that the we of Bpim must be equal to or smaller than Wi, 
although the individual zun may be larger than Wi (see sec. 3). 

The top panel of Figure 11 shows the pixel-pixel beam expansions, -B|(eff) °^ equation 
(5-4). The dotted lines are the results of the individual beams, which are denoted here as 
5L-,. The solid line is the result of the combined beam, which is denoted here as Bl.. Also 
plotted are the B^^^^^^ (eq. [3-6]) and B^^^^^^^ (eq. [3-7]). Here we have used the MAXIMA-1 
scans and a pixel size of 5 x 5 square arcminutes. The bottom panel compares all the above 
Bj to -B|£(jnjd) = 2/[-B2Wsm) + ^T.£ims)\ (^1- [^"6]), with all line styles the same as indicated 
in the top panel. 

We first see that all the -B| have close shapes, with considerable discrepancies only at 
high £ where the amplitude of -B| is small. Second, the bottom panel confirms that the B'^,--^ is 
well constrained by B^^.^s and B^^.s (see eqs. [4-22] and [5-9]), whose fractional difference 
is roughly given by VV| (see definition [3-5]), square of the dashed line in Figure 10. Hence 
according to equation (6-7) and Figure 10, we know that the maximum fractional error in 
the final Cg estimates by taking B'^.^^^. = B^^,^^. for the MAXIMA-1 data will be about 

small error, we still take 

£(c) '^^'■^ -^E£(mid) 



W^/2 fs 5% for £ < 2000. Although this is already a small error, we still take -BLg) = -B^) in 
the MAXIMA-1 data analysis for higher accuracy. The difference between Bjj^.^^ and -B| 
is manifested by the non-zero solid line in the bottom panel of Figure 11. 



-31- 






-0.4 (3) 150 GHz 



-0.4 -0.2 0.2 0.4-0.4 -0.2 0.2 0.4 

degree 



Fig. 9. — Iso-height contours of the beams used in the analysis of the MAXIMA-1 data 
(HOO). Contours correspond to the 90%, 50%, 10%, and 1% amphtude levels. 
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Fig. 10. — Indices of asymmetry of MAXIMA-1 beams (dotted lines) and of the noise- 
weighted combination (solid line). Also plotted is the index of combined asymmetry of all 
MAXIMA-1 beams (dashed line). 
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In addition, we have also verified that for both the individual and the combined beams, 
the approximation B'^,f^^ ^ -Sp£(ms) (^l- [^"9]) i^ accurate within 1% error for i < 2000. 
This means that in general situations one can simply use -B„£(ms) ^s the B'j,^-. to avoid the 
complicated procedure of evaluating the -B^f^g) of equation (5-4) with (5-1). 



Using equation (5-11) with the 3'^,^-. replaced with B ^^i^^-J^^,-. (see eq. [8-9]), we 
tested to what extent our formalism biases the CMB angular power spectrum estimate. We 
simulated a CMB signal 7^ in the time domain. Each time-domain point is allocated the 
pointing coordinates of the MAXIMA- 1 scan and the signal is convolved with the measured 
MAXIMA- 1 beams. In the MAXIMA- 1 scan most pixels are scanned in two different direc- 
tions. We then added time domain noise rit which has the MAXIMA- 1 characteristics: an 
overall white noise, with a 1// behavior at low frequencies due to the receiver response and 
a power law at high frequencies due to the electronic filtering. We call this {dt = 'jt + nt) 

simulation (a). We repeated the procedure to generate simulation (b), in which the CMB 

2 

signal is convolved with a symmetric beam whose power spectrum is identical to B^^t^^y 

Both simulations were then analyzed in exactly the same way, using the procedure described 

in section 2, with the approximation (5-9), (5-11), (8-9), and (8-12). Here we have employed 

the quadratic estimator (Bond et al. 1998) to estimate the power spectra C£(a) and C^(b) for 

simulations (a) and (b), respectively. (The quadratic estimator was implemented by two 

independent codes, one of which is that by Borrill (1999) and the other by the first author, 

and yielded consistent results with less than 0.1% discrepancy.) We then use 

A, = ^^(-^~^^(^) , (9-3) 

where rf(a) is the error bar associated with Q(a), to quantify how much our formalism biases 
the Ci estimates. The entire procedure is repeated six times to yield six independent A^. In 
Figure 12 we plot A^ Vs. d, for the six realizations, the means of these six sets of A^, and the 
standard deviations. As we can see, the means are within 10% of the error bar sizes r^(a) of 
each Q(a). 

With the same scan strategy, pixelization scheme, and noise property, we repeated the 
same test using an extremely elliptic Gaussian beam of 5 by 20 arcminutes in FWHM (the 
one we used previously). We found again that the means of A^ are within 10% of the error 
bar sizes rf(a) for l < 2000. We therefore conclude that our formalism does not bias the Ce 
estimates. 

Finally, we consider the uncertainties in the Ce estimates resulting from uncertainties in 
the measurement of the beam shape, as discussed in section 7. HOO quoted an uncertainty 
of e = ±5% in the measurement of the MAXIMA- 1 beams. The dominant contributors 
to this uncertainty are of the type discussed in section 7 (eq. [7-1]) and contribute to an 
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Fig. 11. — Pixel-pixel beam expansions B'j,^-. of MAXIMA- 1 beams and their noise weighted 
combination (top panel), and a comparison of these results (bottom panel; see text). Also 
plotted are the Bl.., . and El 
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Fig. 12. — Results of simulations testing whether our formalism biases the CMB angular 
power spectrum estimate. The simulations use the MAXIMA-1 scan strategy. Plotted is 
the difference between angular power spectra calculated from simulations that have a CMB 
signal convolved with symmetric and MAXIMA-1 (asymmetric) beams. The difference is 
normalized by the errors of one of the power spectra (see the text). The dotted lines show 
results using different realizations of the CMB signal and the noise. The boxes are the 
averages and error bars are the standard deviations. 
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uncertainty in the Ci estimates that is correlated between different d, bins. Substituting 

f 

p£(ins) 



2 

this value into equation (7-5) and using the B .,^^-. we calculated previously, we obtain the 



estimated uncertainties Ac^ = dCijCi in the Ci estimates. Figure 13 shows the results. As 
one can see, the estimated uncertainties in the C^ estimates are |AcJ < 6%, 17%, and 40% 
for i < 500, 1000, and 1500 respectively. When we include the MAXIMA-1 window-functions 
of i, we find that for the bands used in HOO, the beam size uncertainty causes less than 4% 
and 11% uncertainty in the Ci estimates for i < 410 and 785, respectively. 



10. DISCUSSION AND CONCLUSION 

First we summarize the treatment of asymmetric beams for yielding accurate Ce esti- 
mates in CMB anisotropy experiments: 

1. Based on the measured individual beam pattern Bq, one calculates the 10 A, wq£, using 
equation (3-4), to quantify the level of asymmetry of the beam on different angular 
scales. If taQ^/2 is below the tolerated maximum error for the Ce estimates at the i range 
of interest (see eq. [6-5]), then one takes B^.^^^. = B^^,^^^. as given by equation (6-4) 
(see also eqs. [3-2] and [3-3] for definitions), and goes to step 4, otherwise step 2. The 
resulting errors in the Q estimates by taking Bl^^^. = B"^^,^^. is quantified by equation 
(6-5). Similarly, when combining data from photometers of different beam shapes, one 
first calculates the lOCA We using equation (3-5), and then employs condition (6-7) 
for the same check. If W|/2 is small, then one takes -B|(gfj) = -B|^(mid) ^^ given by 
equation (6-6) and goes to step 4, otherwise step 2. 

2. One checks if all the pixels have the same shape, and if the time-stream beam Bt 
remains unchanged throughout the entire observation. If either or both of these hold, 
then one goes to step 3. Otherwise, one calculates -Bn£(cff) ~ ^npehns) using equations 
(4-11), (4-12), (8-5), and (8-8), and then go to step 5. 

3. One calculates the average pixel-beam expansion Bpim using equations (4-10), (4-11), 
and (4-12). We note that equation (4-10) also works for combining data sets from 
different photometers with different beam shapes, as long as the noise level Ht is well 
taken into account. One then calculates the power spectrum -Bp£(ms) of Bp£m (see eq. [4- 
15]). This can be implemented using the form of equation (4-17) to save computation 
time, i.e. one calculates the weighting function /(/9) first, with discretized P, and then 

2 

the -Bp£(ins) accordingly. A useful check of this result is provided by equation (4-20) or 
(4-21). One thus takes Bj^^^. ^ -Sp£(ms) according to equation (5-9), and goes to the 
next step. 
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4. To incorporate the smoothing effect due to the pixehzation of the map, one employs 
equation (8-10) to obtain the pixel-pixel beam expansion B^^.^^^. ^ ^e(cs)^'e(ms)- ^^ 
general, the associated H^^ ^-^ can be obtained by multipole transforming the n(x) 
that is defined in equation (8-7). If all pixels are regular squares, one can instead use 
the convenient result in equation (8-12). 

5. One then employs equation (2-5) to make a map, and equations (2-6) (or alternatives 
like the quadratic estimator), (2-7), (2-9), and (5-11) to estimate the ^-banded power 
spectrum Cb- We note that in equation (5-11), one replaces the B'^.^g. with the -Bn£(cff) 
obtained previously. 

6. The uncertainties in the final band power Cb resulting from the uncertainties in the 
beam measurement can then be calculated using equation (7-5). In cases where the 
beam has a Gaussian form, one can instead use equation (7-8) with the condition (7- 
10) to estimate the uncertainties. These uncertainties need to be incorporated in both 
the final Ce estimates and the estimates of cosmological parameters. 

In previous sections, we developed the above treatment for asymmetric beams in order 
to obtain accurate Ci estimates at smaller angular scales. This treatment employs the 
symmetric-beam approximation, where the originally asymmetric beams are symmetrized. 
The smoothing effects due to the pixehzation of the CMB map are taken into account. The 
resulting uncertainties in the Ce estimates due to the uncertainties in the beam measurement 
are also estimated. In addition, we derived the conditions under which one needs to employ 
this formalism to account for the asymmetry of beams. We demonstrated certain key points 
by using a simulated highly elliptic beam, and the beams and data of the MAXIMA- 1 
experiment, where the asymmetry is mild. In particular, we showed that in both cases the 
formalism does not bias the final Ce estimates. 

In spite of the power of the new formalism in dealing with various practical situations 
where the beams are not symmetric, we should note that it may break down under certain 
circumstances. First, if the sky patch to be analyzed has an extremely irregular shape, 
then the important result B'^,f^^ ^ -Sp£(ms) (^l- [5-9]) may be invalid due to the nonuniform 
distribution of (/? at each given Ax (see eqs. [5-6] and [5-7]). Nevertheless, the formalism as 
a whole is still valid in this case, because one can instead employ equation (5-4), B'^^^^s = 
{^pp'e)y although it is more computationally expensive. Second, if the total numbers of the 
pixels (A/^) and of the temporal samples {N't) are not large, then some statistical averages 
taken in the formalism may not be appropriate (e.g., eqs. [4-8], [5-4], [5-9], and [8-5]). This 
will cause the violation of some main results like equations (5-9), (8-8), and (8-10). However, 
since the A/^ and A/f are not large in this case, one can always employ the full treatment 
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of asymmetric beams as described by equation (2-8). The main results of our formalism 
are needed only when Mp and Mt are large enough to cause computational difficulty in 
implementing equation (2-8). We note that even in the full treatment of asymmetric beams, 
our results in dealing with the extra convolution effects due to the pixelization of the map (see 
sec. 8) can still be employed. Third, the main results of our formalism have assumed that the 
experimental noise in the temporal samples is independent from each other (i.e., the white- 
noise assumption; see eq. [4-4]), so these results may not be suitable for experiments that 
have strongly correlated noise. Nevertheless, as argued in the appendix, most experiments 
should have only mild departure from the white noise, and this departure does not affect 
our main results. In general, one can use condition (AlO) or equation (All) to choose a 
proper pixel size, so that the white-noise approximation is still appropriate. As we have 
also numerically verified, our formalism does not induce any bias in the final C^ estimates in 
the presence of the nonwhite noise in the MAXIMA- 1 data. Even if the experimental noise 
is extremely nonwhite, we can still deal with asymmetric beams by employing the general 
results in our formalism. This means the use of equation (8-2), together with equations (8-8) 
and (5-11) for the C^ estimation. 

In conclusion, we have proposed a complete and well justified formalism for the data 
analysis of CMB anisotropy experiments. This formalism is very fiexible and therefore well 
suited to a wide spectrum of circumstances, especially when the experimental beams are not 
symmetric. No matter how irregular the beams are, the formalism always provides a both 
computationally economical and statistically plausible way to estimate the angular power 
spectrum of the CMB. We expect this formalism to be useful not only for the small-field 
experiments, but also for the full-sky experiments like PLANCK and MAP. 
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A. Non-white noise 



In this appendix, we consider the pixel-beam expansion Bp£m (eq. [4-3]) in the case 
where the noise is not white (violation of eq. [4-4]), i.e., when it is correlated between pixels. 
We shall show that even in this case, for the purposes of determining Bp^m and thus the 
pixel-pixel beam expansion, the white-noise approximation (4-4) is still appropriate under 
certain conditions, which are generally satisfied by practical situations. 

We start with the general requirement (4-3). In the real space, this equation is equivalent 
to 

Al,,N,7lAtpBp{^) = Al,N^lBt{^). (Al) 

We note that the center of the Bt{x.) here is not at x = 0, but at the location x^ of the 
temporal sample at time t. The similar also applies to the -Bp(x). Thus we see that with 
a given pixelization scheme (and thus given Atp), the relation between Bp and Bt depends 
only on the property of Nft- 

In most experiments, the temporal Fourier transform N{f) of Nft in the frequency / 
domain usually has the following structure: a '1//' behavior below /i due to the receiver 
response, a power law above /h due to the electronic filtering, and a white noise of amplitude 
fi^ (c.f. eq. [4-4]) between /i and /h. Since the temporal Fourier transform A~^(/) of N^^ is 
simply the inverse of N{f), we can approximate a usual N~^{f) as 

iV^H/) ^ f^-' [H{f; /h) - H{f; /O] , (A2) 

where /x^ is the amplitude of the white noise part in N{f), and H{f;fn) {n = h, 1) is a 
top-hat window function: 

H(f:f.) = il. !"!{!-{"■« = h. 1, (A3) 

Thus in the real space we have 



for|/|>/„ 

N^t^N-},-N-,], (A4) 



where 

Nj^ = fi-^ sinc(47r|t - t'|/n), n = h, 1. (A5) 

Here we have used the usual definition sinc(a;) = 2 sin(x/2)/a;. Therefore, to test if the white- 
noise approximation (see eqs. [4-4] and [8-3]) is appropriate, we can substitute N^^if- and N^^,] 
separately as the N^l into equation (Al), and see if the resulting equation is consistent with 
the result (8-3). We may thus derive the conditions under which the result (8-3) is a good 
approximation even if the noise is not white. 



-38- 



We first consider the first term N^^i^ in equation (A4), i.e. the effect from the high- 
frequency cut at /h. In the white-noise case, the first term becomes \im.f^^^oo ^u^t^ which is 
a Dirac Delta along the t and t' directions with a peak centered at t = t'. The fact that the 
width of this peak is zero means that the correlation time of the temporal scan is zero, so that 
each sample is independent and a pixel is related to only the temporal hits inside the pixel, 
giving the form (8-3). On the other hand, in cases where /h is finite, the width of this central 
peak (the distance between the first zeros of N^-t}^ from t = t' along the t or t' direction) is 
broadened (when compared to the white- noise case) from zero to 5th = l//h (see eq. [A5]). 
This means that when the noise is white but with a cut-off beyond /h, the correlation time 
of the temporal samples will be increased from zero to the order of 5th = l//h- Therefore, 
as long as the correlation time 5th = l//h is well below the time required to scan across a 
pixel, the pixel will have no significant correlation with the temporal hits that are outside 
the pixel. In other words, the white-noise approximation (8-3) holds as long as 



6Xr, 5th 

— - > — ^ 
6xt 6t 



/h5t' 



(A6) 



where Sxp is the pixel size, 6xt is the spacing on the sky of the temporal samples, and 6t is 
the integration time of each sample. Although 6xt is not a constant in general, its order in 
a single experiment normally remains the same. 

Following a similar line of logic, we now consider the second term N^,] in equation (A4), 
i.e. the effect from the low-frequency cut at /i. In the white-noise case, the second term 
becomes limjj. 



+oA^it/t = 1^ ^, which is a constant along both the t and t' directions. This 



allows us to simplify equation (Al) as 



^Mep5p(x) 



Eb. 



t X 



(A7) 



where N'tep is the number of temporal samples in the pixel p. Thus we see that equation 
(8-3) automatically fulfills the above requirement. On the other hand, in cases where /i is 
finite, the N^,] will remain constant along the t or t' direction from t = t' out to about 
\t — t'\ = 6t\ = l/2/i (the first zeros), beyond which it begins to decay away as power law 
with oscillations (see eq. [A5]). This means that when /i is not zero but finite, the global 
requirement (A7) will be localized as 



E 

t'ep' 



|t-t'|<<5ti 



E 

t'ep' 



|t-t'|<<5ti 



(AJ 



where Bp^t is the pixel beam of a pixel that covers x^. In the summations over t above (the 
summations inside the brackets), we have ignored the contribution from |t — 1'| > 5ti because 
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the amplitude of N^,] decays as power law and the central (and maximum) amplitudes of 
the beams are always unity (i.e. the contribution from |t — t'| > 6ti decays as a power law; 
see eq. [A5]). Therefore, for equation (A8) to hold for the white-noise result (8-3), we require 
the 6ti to be much larger than the time required to scan through a pixel, i.e. 

5t 2m 5xt' ^ ' 

To sum up, we know that even if the noise is not white, the white-noise approximation 
(8-3) is still appropriate as long as 

In general, the 5xt and 5t are given by experiments, and the 5xp is specified by the pixelization 
scheme. Therefore, since we normally have /h ^ /i in experiments, condition (AlO) can be 
easily satisfied by choosing the right pixel size 5xp. A naive choice will be 

but in fact one would usually like to choose a pixel size closer to the smallest limit to fully 
take advantage of the experimental data. In conclusion, one should choose a pixelization 
scheme whose pixel sizes satisfy condition (AlO) (or have an order given by eq. [All]), so 
that one can employ the white-noise approximation, which leads to some main results of 
this paper (see e.g., eqs. [4-10], [8-5], [8-7], [8-10], and [8-12]). If condition (AlO) can not 
be fulfilled (i.e. when the noise spectrum is far from white or /h ^ /i, which is unlikely to 
be the case), one can still use the general results (8-2) and (8-8) to calculate the pixel-pixel 
beam expansion. 
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Fig. 13. — Estimated uncertainties, A,,^ = rfC^/C^, in the C^ estimates resulting from the 
beam measurement. 



